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Recently, Zou, Janel and Schatz (referred to as ZJS 
below) have described remarkably narrow plasmon reso- 
nances in linear arrays of silver nanospheres [l| . Without 
questioning the novelty and significance of these results, 
I would like to point out that the above-referenced paper 
contains two incorrect statements. 

The first statement is about my previous work. 
Namely, ZJS write that in a previous study [2J I have con- 
sidered "...infinite one-dimensional (ID) arrays in the 
quasistatic approximation" . In fact, there was no qua- 
sistatic approximation made in Ref. |2( • The approxima- 
tion that was made was the dipole approximation. These 
two approximations arc distinctly different. For example, 
even in the electrostatic limit, the dipole approximation 
is grossly inaccurate for two touching conducting spheres 
excited by a constant external electric field parallel to 
the axis connecting the spheres' centers Q . On the other 
hand, electromagnetic interaction of small impurities in 
a crystal or of dye molecules in large molecular aggre- 
gates |4| can not be understood within the quasistatics, 
although the dipole approximation may be very accurate 
in this case. 

Perhaps, the source of confusion is that in Section 2.1 
of Ref. [2] I wrote "The object under investigation is a 
linear infinite chain with step a consisting of point-like 
dipole units (monomers)...". Also, in the introduction 
of Ref. 2] , I have suggested that the physical system to 
which the considered model is applicable is a molecular 
aggregate. Later, in Section 5, I have considered a par- 
ticular example in which the polarizability of a dipole, 
a, was given by the quasistatic polarizability of a small 
sphere with the appropriate radiative correction. How- 
ever, the theoretical formalism of Ref. did not put any 
restrictions on a. And, regardless of the form of a, the 
interaction of dipoles was described with full account of 
retardation effects. 

In fact, ZJS also work in the dipole approximation, al- 
though they validate their results by comparison with a 
more general T-matrix solutions. The situation is some- 
what more complicated, however, because ZJS use, in ad- 
dition, an approximation proposed by Doyle in 1989 5] 
in the context of effective-medium theory of the so-called 
extended Maxwell-Garnett composites, i.e., composites 
in which inclusions are not small compared to the wave- 
length. More specifically, Doyle has studied electromag- 
netic properties of a homogeneous host with randomly 
distributed spherical inclusions. The essence of the ap- 
proximation is to consider only dipolc-dipole interactions 
of the inclusions but to assign them dynamic dipole po- 
larizability a. The latter is given by formula @ below; it 



is defined as the linear coefficient between the amplitude 
of incident plane wave and the total dipole moment of 
polarizablc sphere of arbitrary size (assuming, the sphere 
is isolated) and, in that sense, is exact. It can be seen 
that the Doyle's approach only concerns the choice of a 
within the dipole approximation. Thus, it is fully consis- 
tent with the general formalism developed in Ref. 0. 

It should be noted that the accuracy and limits of ap- 
plicability of the Doyle's approximation have not been 
systematically investigated. In one critical study of 
extended Maxwell-Garnett composites Ruppin has 
shown that the Doyle's approximation is consistent with 
the asymptotes obtained in the limit of small volume frac- 
tion of inclusions, and, in that limit, allows one to con- 
sider inclusions with size parameters of at least x ~ 0.5. 
Thus, the Doyle's approximation can be useful for moder- 
ate size parameters. However, if the spherical inclusions 
are in close proximity of each other, the secondary scat- 
tered waves incident upon each of them are no longer 
plane waves. But the dynamic polarizability used by 
Doyle is exact only with respect to incident plane waves. 
Besides, coupling of higher multipole modes excited in 
spherical inclusions can become significant. Therefore, it 
is quite obvious that the use of Doyle'a approximation 
does not fix, in principle, the deficiencies of the dipole 
approximation. 

The second statement concerns the possibility of 
cancellation of the imaginary part of denominator in the 
expression P = aEo/(l — aS) (Eq. 5, or, in a more 
specific form, Eq. (7) of Ref. jj). This is discussed on 
p. 10874 of Ref. ZJS consider the case when the in- 
cident wave vector is perpendicular to a linear chain of 
polarizable dipoles with the period D. The polarization 
of the incident wave is also perpendicular to the chain. 
It is stated that the resonance width, which is related 
in Ref. 1] to the imaginary part of the denominator of 
the above equation, vanishes when 7 < 8tt 3 A/D 3 , where 
7 and A are parameters which specify the polarizabil- 
ity of an isolated sphere. Namely, ZJS use the formula 
a = —A/(lu — uj p + 17), where uj is frequency of inci- 
dent radiation, oj p is the surface plasmon frequency, and 
7 - the relaxation parameter. Assuming that the result 
ImS* = — 87T 3 /D 3 , which is given in Ref. |l| for A slightly 
larger than the interparticle distance D, is correct, one 
immediately can see that the cancellation takes place ex- 
actly at 7 = 8tt 3 A/D 3 . For smaller values of 7, the imag- 
inary part of the denominator becomes, in fact, negative. 
Such result clearly contradicts conservation of energy and 
is unphysical. It was obtained in Ref. due to several 
mistakes which are discussed below. 
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It is convenient to rewrite Eq. (5) of Ref. pj as 



P = 



1/a-S 



(1) 



Given the specific form a — —Aj (u> — lu p + ry), this ex- 
pression differs from Eq. (7) of Ref. [| only by dividing 
the numerator and denominator by the real constant A. 
The quantity S here is the "dipole sum" - an eigenvalue 
of the electromagnetic state of the dipole chain which is 
excited by incident radiation. The imaginary part of the 
denominator of defines total relaxation. 

Note that Im(l/o!) can contain two contributions 
which correspond to absorptive and radiative relaxation. 
Both are strictly negative. On the other hand, imagi- 
nary part of S has nothing to do with absorptive losses, 
since S does not depend on material properties. Thus, 
ImS can only influence radiative relaxation and can be 
either positive or negative. In the first case, the radia- 
tive relaxation is increased compared to that of an iso- 
lated sphere, while in the latter case it is reduced. It 
is important to note that I/a and S satisfy the follow- 
ing general inequalities: Im(l/a) < — 2fc 3 /3 0,0 an d 
ImS* > — 2fc 3 /3 [3, where k = 2tt/\ is the wavenumber. 
Both inequalities follow from the very general considera- 
tion of energy conservation. At the very least, they show 
that the imaginary part of the denominator of Eq. 
can not become negative. The radiative relaxation is 
canceled if ImS* = — 2fc 3 /3 (this possibility is discussed 
below). If, in addition, Im(l/a) = — 2fc 3 /3, total relax- 
ation is equal to zero. Physically, this can not happen 
due to small absorption which is always present even in 
highly transparent materials, deviations from the dipole 
approximation, etc. 

Let us re- write the above inequalities for Aw D, 
which is the situation considered in Ref. 0]. We obtain 
Im(l/a) < -16tt 3 /3D 3 and ImS > -16vr 3 /3£> 3 . The re- 
sult adduced in Ref. 0, namely, ImS* = — 8n 3 /D 3 m — k 3 , 
clearly contradicts the second inequality. This is due 
to two reasons. First, it is incorrect that the far-field 
term Y^jjti(k 2e ni / r ij) dominates the dipole sum S for 
A w D, as is stated in Ref. This would be only 

true for the real part of S. Second, even if only the far 
field term is used in the calculation of S, the result ad- 
duced in Ref. [J is off by the factor of 2. The correct 
contribution to ImS" which comes from the far-zone term 
is [sgn(D - A)]4tt 3 /D 3 . The contribut ion which comes 
from the intermediate-zone term is 27r 3 /3D 3 . The con- 
tribution from the near-zone term is zero. Thus, we have 
ImS = -I0tt 3 /3D 3 for D < A and ImS* = 147r 3 /3L> 3 
for D > A (all calculations are done for D — A <C D). 
It can be seen that the inequality ImS > — 167r 3 /3Z? 3 
is satisfied strongly. Therefore, not only the imag- 
inary part of the denominator can not become 
negative, but its exact cancellation is also im- 
possible in the considered geometry. The small- 
est possible value of Im(l/a — S) is equal to — 27r 3 /3Z? 3 . 
However, it is correct that the radiative relaxation is 



changed by a significant factor when A — D changes sign. 
Thus, -[167r 3 /3L> 3 + ImS] = -I0it 3 /D 3 for A < D and 
-[167r 3 /3D 3 +ImS'] = -2tt 3 / D 3 for A > D, a drop by the 
factor of 5. This can be practically important if radiative 
losses are dominant over absorptive losses. 

Next, we discuss the inequality Im(l/a) < —2k 3 /3. 
This inequality ensures that the dipole contribution to 
the absorption cross section of a particle is not nega- 
tive. It must hold even for nonabsorbing particles and, 
in particular, for 7 = 0. In the case of a small parti- 
cle, this inequality is satisfied if one uses the quasistatic 
polarizability with the inclusion of the radiative reac- 
tion correction: a — a^ s ^/(I — 2ik 3 a^ s ^ /3). Here 
a is the polarizability with the radiative correction and 
qAQ-S) = i? 3 ( e — l)/(e + 2) is the quasistatic polarizability, 
R being the sphere radius. The importance of the radia- 
tive correction is discussed, for example, in Ref. 0, and 
the authors of Ref. Q are also aware of it (see Ref. 0, 
Eqs. 16-18). The expression a = —A/ (u>—u! p +ij) used in 
Ref. does not contain the radiative correction. There- 
fore, its use (together with an incorrect expression for S) 
leads to unphysical results in the limit 7 — > 0, such as the 
total cancellation of relaxation or negative relaxation. It 
should be also noted that the dynamic expression for a 
which ZJS used in numerical simulations (according to 
the Doyle's approximation) also satisfies the above in- 
equality. Indeed, if we take 



3i rmpiimkRW^kR) - ip^kR^'^mkR) 
2k 3 mi>i(mkR)£i(kR) - Zi(kR)4>[(mkR) 



(2) 



where tpi and £1 are the Riccati-Bessel functions, m — 
y/~e is the complex refractive index of the spheres, then 
the Taylor expansion of Im(l/a) in powers of the wave 
number reads 



Im(l/a) = - 
3/c 4 i?(8 



2fc 3 



3Ime 



3fc 2 Ime 



i? 3 |e-l| 2 5i?|e-I| 2 
- 2Ree)Ime 



0(k 6 R 3 ^ 



350 e - 1 



(3) 



The expansion beyond the third order contains only even 
powers of k and it can be verified that each term in the 
expansion is non-positive. The exact equality Im(l/a) = 
— 2fc 3 /3 takes place only for non-absorbing materials with 
Ime = (which do not occur in nature). 

Finally, we discuss the possibility of exact cancellation 
of the radiative relaxation. Note that (i) only the radia- 
tive part of relaxation can be zero, (ii) the total relaxation 
is always nonzero due to nonzero absorption, but can be- 
come, in principle, arbitrarily small, and (iii) such can- 
cellation can not take place in the geometry considered in 
Ref. . Generally, there can be two reasons for cancella- 
tion of the radiative relaxation. The first is symmetry 0. 
Within the dipole approximation, the cancellation takes 
place when the symmetry of a particular excitation mode 
is such that dipole radiation is forbidden. A non-zero ra- 
diative relaxation can still result from higher-multipole 
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radiation, similarly to non-zero decay rates of excited 
atomic states whose decay is dipole-forbidden. The sec- 
ond reason is when photon emission is prohibited by con- 
servation laws, such as the light cone condition (ll). In 
a linear chain of dipole-polarizable particles the cancel- 



lation of radiative relaxation can take place when the 
incident wave vector is parallel to the chain. However, 
the radiative relaxation is always nonzero for normal in- 
cidence. 
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